The characterization of twenty sequenced human genomes.

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten "case" genomes from individuals with severe hemophilia A and ten "control" genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.

A widespread chromosomal inversion polymorphism contributes to a major life-history transition, local adaptation, and reproductive isolation.

The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements.

Common genetic variation and the control of HIV-1 in humans.

To extend the understanding of host genetic determinants of HIV-1 control, we performed a genome-wide association study in a cohort of 2,554 infected Caucasian subjects. The study was powered to detect common genetic variants explaining down to 1.3% of the variability in viral load at set point. We provide overwhelming confirmation of three associations previously reported in a genome-wide study and show further independent effects of both common and rare variants in the Major Histocompatibility Complex region (MHC). We also examined the polymorphisms reported in previous candidate gene studies and fail to support a role for any variant outside of the MHC or the chemokine receptor cluster on chromosome 3. In addition, we evaluated functional variants, copy-number polymorphisms, epistatic interactions, and biological pathways. This study thus represents a comprehensive assessment of common human genetic variation in HIV-1 control in Caucasians.

Genetic variants in IGF-I, IGF-II, IGFBP-3, and adiponectin genes and colon cancer risk in African Americans and Whites.

PURPOSE: Evaluating genetic susceptibility may clarify effects of known environmental factors and also identify individuals at high risk. We evaluated the association of four insulin-related pathway gene polymorphisms in insulin-like growth factor-1 (IGF-I) (CA)( n ) repeat, insulin-like growth factor-2 (IGF-II) (rs680), insulin-like growth factor-binding protein-3 (IGFBP-3) (rs2854744), and adiponectin (APM1 rs1501299) with colon cancer risk, as well as relationships with circulating IGF-I, IGF-II, IGFBP-3, and C-peptide in a population-based study. METHODS: Participants were African Americans (231 cases and 306 controls) and Whites (297 cases, 530 controls). Consenting subjects provided blood specimens and lifestyle/diet information. Genotyping for all genes except IGF-I was performed by the 5'-exonuclease (Taqman) assay. The IGF-I (CA)(n) repeat was assayed by PCR and fragment analysis. Circulating proteins were measured by enzyme immunoassays. Odds ratios (ORs) and 95 % confidence intervals (CIs) were calculated by logistic regression. RESULTS: The IGF-I (CA)( 19 ) repeat was higher in White controls (50 %) than African American controls (31 %). Whites homozygous for the IGF-I (CA)(19) repeat had a nearly twofold increase in risk of colon cancer (OR = 1.77; 95 % CI = 1.15-2.73), but not African Americans (OR = 0.73, 95 % CI 0.50-1.51). We observed an inverse association between the IGF-II Apa1 A-variant and colon cancer risk (OR = 0.49, 95 % CI 0.28-0.88) in Whites only. Carrying the IGFBP-3 variant alleles was associated with lower IGFBP-3 protein levels, a difference most pronounced in Whites (p-trend <0.05). CONCLUSIONS: These results support an association between insulin pathway-related genes and elevated colon cancer risk in Whites but not in African Americans.

A functional polymorphism in the 5HTR2C gene associated with stress responses also predicts incident cardiovascular events.

Previously we have shown that a functional nonsynonymous single nucleotide polymorphism (rs6318) of the 5HTR2C gene located on the X-chromosome is associated with hypothalamic-pituitary-adrenal axis response to a stress recall task, and with endophenotypes associated with cardiovascular disease (CVD). These findings suggest that individuals carrying the rs6318 Ser23 C allele will be at higher risk for CVD compared to Cys23 G allele carriers. The present study examined allelic variation in rs6318 as a predictor of coronary artery disease (CAD) severity and a composite endpoint of all-cause mortality or myocardial infarction (MI) among Caucasian participants consecutively recruited through the cardiac catheterization laboratory at Duke University Hospital (Durham, NC) as part of the CATHGEN biorepository. Study population consisted of 6,126 Caucasian participants (4,036 [65.9%] males and 2,090 [34.1%] females). A total of 1,769 events occurred (1,544 deaths and 225 MIs; median follow-up time = 5.3 years, interquartile range = 3.3-8.2). Unadjusted Cox time-to-event regression models showed, compared to Cys23 G carriers, males hemizygous for Ser23 C and females homozygous for Ser23C were at increased risk for the composite endpoint of all-cause death or MI: Hazard Ratio (HR) = 1.47, 95% confidence interval (CI) = 1.17, 1.84, p = .0008. Adjusting for age, rs6318 genotype was not related to body mass index, diabetes, hypertension, dyslipidemia, smoking history, number of diseased coronary arteries, or left ventricular ejection fraction in either males or females. After adjustment for these covariates the estimate for the two Ser23 C groups was modestly attenuated, but remained statistically significant: HR = 1.38, 95% CI = 1.10, 1.73, p = .005. These findings suggest that this functional polymorphism of the 5HTR2C gene is associated with increased risk for CVD mortality and morbidity, but this association is apparently not explained by the association of rs6318 with traditional risk factors or conventional markers of atherosclerotic disease.

Can Genetics Predict Response to Complex Behavioral Interventions? Evidence from a Genetic Analysis of the Fast Track Randomized Control Trial.

Early interventions are a preferred method for addressing behavioral problems in high-risk children, but often have only modest effects. Identifying sources of variation in intervention effects can suggest means to improve efficiency. One potential source of such variation is the genome. We conducted a genetic analysis of the Fast Track randomized control trial, a 10-year-long intervention to prevent high-risk kindergarteners from developing adult externalizing problems including substance abuse and antisocial behavior. We tested whether variants of the glucocorticoid receptor gene NR3C1 were associated with differences in response to the Fast Track intervention. We found that in European-American children, a variant of NR3C1 identified by the single-nucleotide polymorphism rs10482672 was associated with increased risk for externalizing psychopathology in control group children and decreased risk for externalizing psychopathology in intervention group children. Variation in NR3C1 measured in this study was not associated with differential intervention response in African-American children. We discuss implications for efforts to prevent externalizing problems in high-risk children and for public policy in the genomic era.

Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays.

Determination of copy number variants (CNVs) inferred in genome wide single nucleotide polymorphism arrays has shown increasing utility in genetic variant disease associations. Several CNV detection methods are available, but differences in CNV call thresholds and characteristics exist. We evaluated the relative performance of seven methods: circular binary segmentation, CNVFinder, cnvPartition, gain and loss of DNA, Nexus algorithms, PennCNV and QuantiSNP. Tested data included real and simulated Illumina HumHap 550 data from the Singapore cohort study of the risk factors for Myopia (SCORM) and simulated data from Affymetrix 6.0 and platform-independent distributions. The normalized singleton ratio (NSR) is proposed as a metric for parameter optimization before enacting full analysis. We used 10 SCORM samples for optimizing parameter settings for each method and then evaluated method performance at optimal parameters using 100 SCORM samples. The statistical power, false positive rates, and receiver operating characteristic (ROC) curve residuals were evaluated by simulation studies. Optimal parameters, as determined by NSR and ROC curve residuals, were consistent across datasets. QuantiSNP outperformed other methods based on ROC curve residuals over most datasets. Nexus Rank and SNPRank have low specificity and high power. Nexus Rank calls oversized CNVs. PennCNV detects one of the fewest numbers of CNVs.

Genetic variants in the TEP1 gene are associated with prostate cancer risk and recurrence.

BACKGROUND: Telomere-related genes play an important role in carcinogenesis and progression of prostate cancer (PCa). It is not fully understood whether genetic variations in telomere-related genes are associated with development and progression in PCa patients. METHODS: Six potentially functional single-nucleotide polymorphisms (SNPs) of three key telomere-related genes were evaluated in 1015 PCa cases and 1052 cancer-free controls, to test their associations with risk of PCa. Among 426 PCa patients who underwent radical prostatectomy (RP), the prognostic significance of the studied SNPs on biochemical recurrence (BCR) was also assessed using the Kaplan-Meier analysis and Cox proportional hazards regression model. The relative telomere lengths (RTLs) were measured in peripheral blood leukocytes using real-time PCR in the RP patients. RESULTS: TEP1 rs1760904 AG/AA genotypes were significantly associated with a decreased risk of PCa (odds ratio (OR): 0.77, 95% confidence interval (CI): 0.64-0.93, P=0.005) compared with the GG genotype. By using median RTL as a cutoff level, RP patients with TEP1 rs1760904 AG/AA genotypes tended to have a longer RTL than those with the GG genotype (OR: 1.55, 95% CI: 1.04-2.30, P=0.031). A significant interaction between TEP1 rs1713418 and age in modifying PCa risk was observed (P=0.005). After adjustment for clinicopathologic risk factors, the presence of heterozygotes or rare homozygotes of TEP1 rs1760904 and TNKS2 rs1539042 were associated with BCR in the RP cohorts (hazard ratio: 0.53, 95% CI: 0.36-0.79, P=0.002 and hazard ratio: 1.67, 95% CI: 1.07-2.48, P=0.017, respectively). CONCLUSIONS: These data suggest that genetic variations in the TEP1 gene may be biomarkers for risk of PCa and BCR after RP.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Association of Common Genetic Polymorphisms with Melanoma Patient IL-12p40 Blood Levels, Risk, and Outcomes.

Recent investigation has identified association of IL-12p40 blood levels with melanoma recurrence and patient survival. No studies have investigated associations of single-nucleotide polymorphisms (SNPs) with melanoma patient IL-12p40 blood levels or their potential contributions to melanoma susceptibility or patient outcome. In the current study, 818,237 SNPs were available for 1,804 melanoma cases and 1,026 controls. IL-12p40 blood levels were assessed among 573 cases (discovery), 249 cases (case validation), and 299 controls (control validation). SNPs were evaluated for association with log[IL-12p40] levels in the discovery data set and replicated in two validation data sets, and significant SNPs were assessed for association with melanoma susceptibility and patient outcomes. The most significant SNP associated with log[IL-12p40] was in the IL-12B gene region (rs6897260, combined P=9.26 × 10(-38)); this single variant explained 13.1% of variability in log[IL-12p40]. The most significant SNP in EBF1 was rs6895454 (combined P=2.24 × 10(-9)). A marker in IL12B was associated with melanoma susceptibility (rs3213119, multivariate P=0.0499; OR=1.50, 95% CI 1.00-2.24), whereas a marker in EBF1 was associated with melanoma-specific survival in advanced-stage patients (rs10515789, multivariate P=0.02; HR=1.93, 95% CI 1.11-3.35). Both EBF1 and IL12B strongly regulate IL-12p40 blood levels, and IL-12p40 polymorphisms may contribute to melanoma susceptibility and influence patient outcome.

Genetic Variants of the MDM2 Gene Are Predictive of Treatment-Related Toxicities and Overall Survival in Patients With Advanced NSCLC.

INTRODUCTION: Platinum agents can cause the formation of DNA adducts and induce apoptosis to eliminate tumor cells. The aim of the present study was to investigate the influence of genetic variants of MDM2 on chemotherapy-related toxicities and clinical outcomes in patients with advanced non-small-cell lung cancer (NSCLC). MATERIALS AND METHODS: We recruited 663 patients with advanced NSCLC who had been treated with first-line platinum-based chemotherapy. Five tagging single nucleotide polymorphisms (SNPs) in MDM2 were genotyped in these patients. The associations of these SNPs with clinical toxicities and outcomes were evaluated using logistic regression and Cox regression analyses. RESULTS: Two SNPs (rs1470383 and rs1690924) showed significant associations with chemotherapy-related toxicities (ie, overall, hematologic, and gastrointestinal toxicity). Compared with the wild genotype AA carriers, patients with the GG genotype of rs1470383 had an increased risk of overall toxicity (odds ratio [OR], 3.28; 95% confidence interval [CI], 1.34-8.02; P = .009) and hematologic toxicity (OR, 4.10; 95% CI, 1.73-9.71; P = .001). Likewise, patients with the AG genotype of rs1690924 showed more sensitivity to gastrointestinal toxicity than did those with the wild-type homozygote GG (OR, 2.32; 95% CI, 1.30-4.14; P = .004). Stratified survival analysis revealed significant associations between rs1470383 genotypes and overall survival in patients without overall or hematologic toxicity (P = .007 and P = .0009, respectively). CONCLUSION: The results of our study suggest that SNPs in MDM2 might be used to predict the toxicities of platinum-based chemotherapy and overall survival in patients with advanced NSCLC. Additional validations of the association are warranted.